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^ ' Abstract 

> : 

. The one-loop polarization operator of neutral gluons in the background constant 

I Abelian isotopic, H^, and hypercharge, Hg, chromomagnetic fields combined with 

' Aq electrostatic potential at high temperature is calculated. The case when Aq = 

Cp . is investigated separately. The proper time method is applied. It is found that 

neutral gluons do not acquire magnetic masses in the background fields, in contrast 
to the charged ones. The application of the results are discussed. 

§^ 1 Introduction 

^ ! Investigation of the deconfinement phase of QCD remains of considerable interest for 
: high-energy physics and cosmology. Among the most important objects here is a gluon 



O 



^ ' polarization tensor (PT) containing information on the excitation spectrum of quark- 

gluon plasma. First the QCD PT was calculated and investigated in one-loop order of 
perturbation theory at T 7^ by Kalashnikov and Klimov pp, [2] (see also surveys |SI 
and jl], where the results on higher order contributions are discussed). As it has 
been shown, the space components of the one-loop gluon propagator calculated within 
a standard perturbation theory possesses a fictitious infrared pole at /C4 = 0, A; ~ g'^T 
which could not be removed by any further resummations. These infrared divergencies of 
the termal Green functions provide the most challenging difficalties in understanding the 
internal structure of perturbative finite temperature QCD. It is believed, however, that 
formation of some condensate fields, such as a uniform "colour" magnetic field {He = 
const.) or electrostatic potential (so-called Ao condensate), can improve the infrared 
properties of the theory. These condensate fields may arise in the deconfinement phase of 
QCD due to a peculiar dynamics of non- Abelian gauge fields, as it was argued by several 
authors |ni-[lll- In the paper by Kalashnikov jT^I it was demonstrated, in particular, 

*Permanent address: Dniepropetrovsk National University, St. Naukova 13, Dniepropetrovsk 49050, 
Ukraine. 



that the Aq condensate shifts the fictitious pole and introduces the gluon magntic mass 
of the order ~ g'^T'^. At the same time, in Ref. ^1 it was discovered that in the 
presence of the external Abelian chromomagnetic fields H the transversal charged gluons 
acqure a magnetic mass m^^g^ ~ g'^ VJH T which is generated within the one-loop 
polarization operator. It acts to stabilize the external field. In Ref. it was found 
within the SU{3) gluodynamics that at high temperature a specific combination of the 
Abelian hypercharge, H^, and isotopic spin, H3 , fields is generated and is stable due to 
this magnetic mass. It is also of the order ~ g^T^- The tachyonic (unstable) modes of the 
transversal charged gluons, which appear in the energy spectrum of the charged vector 
particles when the homogenious magnetic field is applied to the system, are removed 
by these high-temperature radiative corrections. Moreover, an imagenary part of the 
effective potantial (EP) of the background fields is cancelled if the contribution of the 
daisy diagrams with this magnetic mass is taken into consideration. Hence, one has to 
believe that the non-trivial configuration of the classical magnetic fields and ifg is 
generated in the deconfinement phase. 

It is interesting to see in actual calculations whether or not the magnetic mass of the 
neutral gluons is generated in the external field at high temperature. Actually, this is 
not expected because on general theoretical grounds the fields belonging to the Abelian 
projection of the non-Abelian groups remain massless. It is also important to know 
whether or not the fictitious pole of the neutral gluons is preserved when a magnetic field 
and Aq is present in the system. 

The aim of the present paper is to calculate the one-loop polarization operator of 
the neutral gluons in SU(3) gluodynamics in the external fields H^, and and Aq{ or 
A4) electrostatic potential at high temperature and check whether the full propagator of 
neutral gluons Q3 and Qs contains the fictitious pole leading to the infrared instability. 
If this is not the case, one is able to conclude that the formation of the condensate fields 
play the role of an infrared regulator and the transversal components of neutral gluons are 
unscreened. It is necessary to note that at zero temperature this problem was investigated 
in Ref. [201 • ^i^^ begin with the case when both the chromomagnetic fields and the 
electrostatic potentials are present in the system. Then the case of gA^ = will be 
separately analysed. We will restrict our consideration to the one-loop approximation. 
To evaluate integrals over a three-dimantion momentum the Fock-Schwinger proper time 
method will be applied. The most essential steps of calculation are given in the Appendix 
1. 



2 Calculation of the polarization operator 

We start our analysis with the expression of the Lagrangian of neutral gluons in (Eu- 
chdean) SU{3) gluodynamics: 

Lneut.,L = " IQ^QU " Ud,Ql){d.Ql) " \{d,Ql){d^Ql) 

+^\l\9iQl + T^Qimi{d,w^, - d^w^^) - {h.c.)) 
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-l9'r,.xpQ',Ql{w,\w,- + W,\W: 



gh- 



Here the following basis of charged gluons Q't (a = 1, 2, 4, 5, 6, 7) 



^{Ql±^Ql)■ (2) 



is introduced. The external potential is chosen in the form B*^ = d'^^B-^^ + 6°'^Bsfj_, where 
i?3^ = H^5fj,2Xi + 5^igA'i and Sg/i = -f^s^^s^i + '^m49'^8- In these formalae the notations 
and correspond accordingly to Aq~^ and Aq""^ electrostatic potentials. The constant 
chromomagnetic fields are chosen to be directed along the third axis of the Euclidean 
space and a = 3 and a = 8 of the colour >S?7c(3)-space: F^/^* = 5''^F^=^ + 5"^F;=^ 
F^2 = "-^21 — o = 3,8. From the Lagrangian one can easely derive the diagrams 
describing propagation of the neutral gluons in the background fields. 

In the one-loop approximation the PO of neural gluons is determined by the standard 
set of diagrams in Fig. 1, where double wavy lines represent the Green function Gr fiu{x, y) 
for the charged gluons, dashed double lines represent the Green function D{x,y) for the 
charged ghost fields. Thin wavy line corresponds to the neutral gluon fields Q^'^. In 
the operator form the above Green functions are given by the expressions (in Feynman's 
gauge) 



1 



D{P) 



1 

p2 ' 



A 



1 ± 



a/6 ^8 



To calculate the PO we make use of the proper time representation and the Schwinger 
operator formalism T7J. The PO of the neutral gluons in the background fields at T 7^ 
can be written as 



where 



n 



a=3 



2^ (27r)3 



n^.(A:,F), 



(3) 
(4) 



n^.(A:,F) 
n;..(A:,F) 



E n;.(A;,p), 



r=2,3 



K-{k, P) = [Vr pc^AP^ k) Gr /3a(P) P, A:) G, aa{P - A;) - 26^,Gr aaiP) + 

(2P - k)^D{P){2P - k),D{P -k)- 25^,D{P)] }, 

r^a,/3 = {2P - k)f,5al3 -2{kaSi3^ - kfjdat^), 
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P4 = 27r/r + gAs, Pi = id, + gB^i for r = 1 and P4 = Svrff + J'^fi±gA^ 



1 nT^r.- fnr = 9 5^ r#^c;r~»<3r'f i^r#^Kr- 7 = -1-1 -1-9 

1 ^3 



Pj = idi + J^X±gB^i for r = 2, 3, respectively; Z = 0, ±1, ±2, . . ., and 



Assume now that the values of potantials A^^ and satisfy the following conditions: 
gA'i <^ T and gA^i, -C T. This is natural because the quantities gA^^^^ are expected to be 
of order g'^T, as it is pointed out in Refs. [ID], [HI for SU (2) case. To investigate the high 
temperature limit of (jH)) and one can take the / = term only in the sum over P4 [Hj. 

To evaluate the expression for the PO let us apply the Schwinger proper-time method 
modified for the case of high temperature. From technical point of veiw, this case is similar 
to the zero temperature one, so one may consult for more details, for example, to Refs. jl8j- 
|22j . where the PO of photon as well as neutral gluon in the external (chromo)magnetic 
field were calculated at T = . The basic steps of the calculating procedure are noted 
in the Appendix. For simplisity it is convenient to introduce the following notations: 

H^ = ^-X^Hs; A^ = ^-^,^As; ^ = ^> = 
Then the final result of evaluation and reads: 



n-4==^'' = -n4T''' = ^m^^%,S, Kf^' = nl'is.S - ^a=3,8. (7) 

Here the quantities n^*l3 8, i = I, ..,4, and ^^0=3,8 are 

niia = n« + iniig, ^a=3 = ^ + l^a=S, (8) 




n 



du 



sh{x] 



+ 4:sh{x) 



(11) 



Sg"^ r°° dx 



sh{x) 



+ 4s/i(x)]{v^e-^'"+ + V^e-^™-}, (12) 
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where 



^ = xu{l - v)^- + k\—^, xu{l-u)- ^ 



and 

_ ch{x) - ch{x{l - 2u)) , 2 _ , 2 , 1,2 

C, — i/\ , — Ki -\- ^2 



Exact expressions for the functions Z*^*-* and /^^ are adduced in the Appendix 2. The 
matrix Bij is a usual two dimension antisymmetric tensor, 

Bij = eij = Si2Sij — 6ii52j- 

The spatial part of the PO is transversal manifestly, as it is required by gauge invariance. 
Note that U^^=^'^) = o for = /l± = 0. 

Now let us consider the high-temperature expansion, gH-s s ^ T"^, (fi'^s.s)^ ^ T"^, of 
the expressions in Eqs.flHj)- (fT^ . Assuming that the quantities gH^^s and {gA^^sY are of 
the same order of magnitude, we investigate the two separate regimes: | ^ g^T and 
I ^ |> gT. In the former case, with the additional condition k'j_ <C ^'-^3,8 and k^ < gH^^g, 
the main contributions to integrals come from the integration domain where x ^ 1. 
Carrying out integrations we obtain 

nii3 = n«(^//3;^;m) + lnli3, (13) 

ntt''^ = l{u^'='^'\gH^; u^; m+) + Il^'='''\gH_; m„)}, (14) 

n28 = ^{^n(3n^i^+;^+;"^+) + ^n(^H^7i^-;^-;rn_)}, (is) 

n^s = l{i^rU^'\9H^; m^) + i^rU^'\gH^; u.; m_)}. (16) 



expressions: 



Here v = -rjr, = . \ and functions n'^*)(Q;; /3; 7) are represented by the following 



^ '^'^^ 27rv/^4(7/5+l) + /?2^ ^ VlT^ ^ 



87rv/^^(7- 1)(7 + /?- 1) 
1+/3 1 



i + /3(7 + /3 + i)^v/y^ VYTWTi 



)]-n«, 



TCy^ [4(7 - 1) + /9] ^ 
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where according to (fT^ - (fT^ instead of variables a, (3 and 7 one has to substitute 
gH^, ly, m or gHh±, m±, respectively. For \E'a=3 and ^1^0=8 we have: 

^'a=3 = "^{gHs; m) + l'^a=8, (17) 



^a=s = ^(*(^^+; m+) + ^{gH_; m_)), (18) 

where 

7) = ^v^[^^ + \/7-i]- 

For the values m = 1 and/or m± = 1 the functions n^*"^'^-' and 11'^^^ become divergent 
whereas n^^-' is equal to zero. 

In the case of | k \> gT and k\ ^ {gA^QY (but m > 1 and 'm± > 1), the main 
contributions to integrals come from the region x ~ 0. Expanding the integrand functions 
into the power series over the veriable x, one can obtain for the spatial components (jH}: 

n^'' ~ (^^. - ^)^'ni23,8 + {BkUBk)j^^:i3,s- (19) 

Here 

ni'i = -|a n23 = 5c. nii = -|La ni2. = ^c. c = ^. (20) 

It is remarkable that the quantities (j^n|) . which are, of cause, only the leading terms of 
perturbative expansion, do not depend upon the condensate fields. For the momentum 
scale A;_|_ ~ T the constant C is of order g"^ and, therefore, perturbative theory is actually 
governed by the parameter g"^. However, for the scale k± ~ gT ^ T the effective 
expantion parmeter becomes g. Hence one can see that perturbative features of the 
model are aggravated with the decreasing 



3 Polarization operator in the extarnal magnetic fields 

In this section we consider the PO in the external chromomagnetic fields H^ ^ (but gA^ g = 
0). We merely put the parameters A3, A±, in the Eqs. fPl- IT^ equal to zero. In this case 
the integrands in the r.h.s. of Eqs. (Pj)- (jl2|) are nonanalytical for large x. To ensure the 
convergence of integrals with respect to x one has to rotate the integration contour by 
the standard rule: x — >■ ix. Then, assuming again that kj_ <C gH^^s and A;| ^ gH3,8, the 
main contributions come from large x and we obtain: 

nli3 = n«(^?^3;z^) + ^n«8, (21) 
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where 



n(^)(«-/?) = -^[ ^ 2(1-2^ 

^ '^^ 27rv/^^4v/:5TT(l + /3(/? + l)) 1 + 4/9 

4(1-/9)^1 + 4/? 4(l + /?(/3 + l))^^ 
The Debye masses of neutral gluons are 

,2 



and 



Re{Ulf) = ^-T JgH, + J gH,-iJ K + V^-) (23) 



4 



Re{W=^) = ^Ty^(^+ ^). (24) 

Quantities (PT|) and (j221), as well as 1144, include imaginary parts reflecting the existence 
of the tachyonic mode in the tree-level spectrum of charged gluons. It should be noted, 
that expressions for 1144 represent the next-to-leading terms. To calculate the leading 
terms one has to perform summation over the discrete frequencies P4. 

4 Discussion 

To discuss the results obtained, let us consider the full propagator of the neutral gluons 
Qa=3,8_ one- loop order the transversal part of the propagator spatial components has 
the following structure 



P [p(i + n(i)) + A;in(2)](i + n«)' ^ ' 

where the functions H^^'^) are given by Eqs.dTHI), (HH), (j2IIl), dill) and (j221)- In the case 
of 5(^3 8 7^ (see Eqs. (fTSj) . (fT^ and ), the full propagator of does not contain 
a non-trivial pole. Hence, one has to conclude that the neutral gluons do not acquire 
magnetic masses in the presence of the background fields A3 § and i^3,8- 

Here a more serious problem arises. Namely, if the condensate fields are of the order 
g'^T, as it was argured in Refs. P-fH], then, for the case of k\ <C JgH^ g (see flTS jl - (fT^ ). 



the factors / appearing in (fT ^ - (IT^ turn out to be of order 0(1) and the perturbative 

expansion breaks down for the momentum scale k± -C g^T. Therefore one cannot explore 
the infrared region {k 0) by usual perturbative methods and our conclusion is valid for 
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the scale kj_ > gT, only. In this region perturbation theory is reliable (see Eqs. (jl9p - (j20j) 
and the text below). 

In the case of the chromomagnetic fields been taken into consideration the quantities 
were found to be complex, and Eq. ()25p can be rewritten as 

_ kikj 1 + i?en(i) - i /mn(y 



(Bk)(Bk) (i + ^enw-^^"^nw)n(^) 

^ [(P(l + i?enW) + A;in(2))2 + (p/mn(i) + A;i/mn(2))2] ' ^^^^ 

This expression has also a pole at A;2 = 0, only. However, the imaginary part that arises 
in Eq. has a "tachyonic" origin, as it was mentioned above. Really, the calculation 
of Iljj (as well as 1144) has been carried out with the bare propagators of the charged 
gluons substituted into internal lines of diagrams. This results in a non-analyticity of 
integrands with respect to the variable x in the 11^,^. In this sence the carried out one- 
loop calculation of the PO appears to be insufficient: to obtain a correct independent 
of the imaginary part expressions for (j26p . the charged gluon propagators accounting for 
the magnetic mass derived in the paper 12j must be used. But now, when we know the 
origin of the imaginary part, it does not matter when the problem on the magnetic mass 
of the neutral gluons is investigated. 

In the present paper it was straightforwardly demonstrated that the transversal neu- 
tral gluon fields are not screened by thermal fiuctuations if the non-trivial condensates 
present in the QCD deconfinement phase. We arrived at the following picture when the 
assumed formation of condensate fields gA^^g and gH^^g determine the effective masses of 
the charged gluons while the neutral spatial components do not acquire magnetic masses 
in the fields. It is resonable to suppose that this picture will be also valid when only 
chromomagnetic fields gH^ g are generated in the system although higher-order contri- 
butions to the neutral gluon PT must be taken into account in this case. It is worth to 
emphasize that in the infrared region, A; — > 0, the full propagator ()25|) does not contain 
the "fictitious" pole. This is in contrast to the case of trivial vacuum [I], 



5 Acknowledgments 

One of us (A.S.) thanks M. Bordag for helpfuU discussions and the graduate college 
Quantenfieldtheorie at the University of Leipzig for support and friendly enviroment. 
This work (V.S.) was supported in part by the grant from DFG No: UKR 427/17/03. 



8 



6 Appendix 1 



To illustrate the basic stages of evaluating the PO (jnHH) let us consider the integral: 



2 . ^3p 



I = 

i5 J (27r)3 



n,,(^,p), 



which represents the contribution of the charged fields (and the corresponding 

ghosts) to the H^j"^ at high temperature. The rest components of the tensor J^,^ are 
calculated analogously. Following the standard procedure we introduce a proper time for 
each propogator appearing in Hijik, P): 



D{P) 




oo 



Jo 

Then, the whole expression for I can be rewritten in the form: 

,2 . ^3p 



kj = ^ dsids2 J j—:^exp[-{si + S2){gA3f] 

{ [Fi a,miP^ k) A„„,(ai) Ti ,,,„(P', k) A,i{cr[) - 2{2P - k),{2P' - k)^] 6^=1} 



d^P r 



TrA{a'l) - 2 



where Ti u^rn = (2P - k)i5im - 2{ki6mi - kmSu) is the vertex factor, 

er=i = e-''^' e-'^^^-~^^\ Aij{x) = Rij - Blch{x) - iBijsh{x) 

and variables ai, a[, a'( are 

ar=i = 2igHsSi, a^^^ = 2igH3S2, a"=i = 2igHss. 

We introduced the following designation: P/ = {exp[—2igF2Si]P)i, Pi = idi + gB^i. The 
matrixes Bij and Bf^ are: 



P 



tj — "jswsj) Bij — 6i2Sij — Sii62j, B-j — BijBij. 



Next, three-dimensional integration with respect to P in Jjj is carried out by means of 
the transition to the conjugate variable X-: 

[X,,P,]=^S,,. 

By using the eigenstates of the operator Xi as determined by the condition X- = 0, the 
integral over P can be represented as 



d^P 
(2^ 



f{P) = {X' = 0\f{P)\X' = 0). 



Hence, performing the following transformation of variables Si and 82'- Si = s{l — u) 
S2 = su, we have for lij : 



9 



1 oo 







9 /-oo 



TrAK)-2expH(^A3)2](e-^l, 



X' = 0) = (...). Now one need 



For convenience we use a notation {X' = | .. 
to calculate the quantities (6*^=1), {PiOr=i) and {PiPj9r=i) according to the procedure 
described in Ref. jTTj. The result reads 



{Pi9r=l)= (^.=1) 



D Ji\D 



k] -ig 

3 

9H3S 



F 



(47rs)3/2 shigH^s) 



$3 



^3 = klsu{l - u) + k"^ ^ 



2gH3' 

where A = e-^'^uFsu -i^ d = e-^'^uFs _ 1^ ^ ^2 ^ ^2 

ch{gHss) - ch{gHs{l - 2u)s) 



c 



Note that 



shigHss) 
1 gHss 



(47rs)3/2 shigHss) ' 
Finally, after integrating by parts, we arrive at: 

1 00 



9' T 



87r3/2 



du 




dx 
sh{x) 



X exp[—^ — xm]Mij{x, u). 



Here designations (jH} are used and x = gH^s. The matrix Mij is defined by 

Mij = {2ch{2x)[{pk)i {pk)j - Pij{kpk) + (AA;)i(Afc),] + S{Bk),{Bh)jCsh{2x) 

+4:[-{^{a')kU^{a)k), - {A{-a)kUA{-a')k)i 
+A,,{a')CkA{a)k) + A^Ji-a)CkA{a')k)]}, 

where a = 2x{l - u), a' = 2xu, p = {I - 2u)R - X = (B,^= ^^^^f|^(^. It can be 
easely verified that quantity lij is manifestly transversal, kilij = kjlij = 0, as it should 
be due to the gauge invariance. 
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Now we can apply described above procedure to evoluate the 11^^ 
j;iven by: 



a=3,l 



The result is 



n 



(0 



du 



^ g— '3?+— a;m+ _j_ ^ ^— — xm- 



VgH- 



u 



values, K/^i — QQ, as 



It is conveniant to rewrite the operators II^,- ' , using their eigenvectors, &f, and eigen- 

5 '^a=3,85 



n: 



a=3,8 _ 

p=l 



where = (Bk)i, = (Rk)i + ^{B'^k)i and = hi. The eigenvectors 6f satisfy 
to the condition of completness: 

V = 5 - 



Hence, since ^ = Q because of transversality of the H"- we obtain (jH]) 



7 Appendix 2 



The functions and l^\x,u) are: 



1 



/(^) = 4[2c/i(x)c/i(x[l - 2m]) - -(1 - 2u)^ch{2x) 



= 4[c/.(2a;) - 2sh{2x)C - ^c/i(2x)(e' - C')] - 



= AgAs \2shi2x) - c/i(2x)cl , = 8igAsfchi2x) 



and 



,(i=l,2) _ r{i=l,2) ,(3) _ ^ .(3) , 



(4) _ (^±)' 



,/(4) 



^ _ - 2u)) ^ _ ch{x) - ch{x{l - 2u)) 



sh{x) 



sh{x) 
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Figure 1: Polarization operator of neutral gluons in the one-loop approximation. 
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